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ABSTRACT 

Context. Stellar winds observed in asymptotic giant branch (AGB) stars are usually attributed to a combination of stellar pulsations 
and radiation pressure on dust. Shock waves triggered by pulsations propagate through the atmosphere, compressing the gas and lifting 
it to cooler regions which creates favourable conditions for grain growth. If sufficient radiative acceleration is exerted on the newly 
formed grains through absorption or scattering of stellar photons, an outflow can be triggered. Strong candidates for wind-driving dust 
species in M-type AGB stars are magnesium silicates (Mg 2 SiC >4 and MgSiCfs). Such grains can form close to the stellar surface, they 
consist of abundant materials and, if they grow to sizes comparable to the wavelength of the stellar flux maximum, they experience 
strong acceleration by photon scattering. 

Aims. The purpose of this study is to investigate if photon scattering on Mg 2 SiC >4 grains can produce realistic outflows for a wide 
range of stellar parameters in M-type AGB stars. 

Methods. We use a frequency-dependent radiation-hydrodynamics code with a detailed description for the growth of Mg 2 SiC >4 grains 
to calculate the first extensive set of time-dependent wind models for M-type AGB stars. This set includes 139 solar-mass models, 
with three different luminosities (5000 L Q , 7000 L G , and 10000L o ) and effective temperatures ranging from 2600 K to 3200 K. The 
resulting wind properties, visual and near-IR photometry and mid-IR spectra are compared with observations. 

Results. We show that the models can produce outflows for a wide range of stellar parameters. We also demonstrate that they 
reproduce observed mass-loss rates and wind velocities, as well as visual and near-IR photometry. However, the current models do 
not show the characteristic silicate features at 10 and 18 pm as a result of the cool temperature of Mg 2 SiC >4 grains in the wind. 
Including a small amount of Fe in the grains further out in the circumstellar envelope will increase the grain temperature and result in 
pronounced silicate features, without significantly affecting the photometry in the visual and near-IR wavelength regions. 
Conclusions. Outflows driven by photon scattering on Mg 2 Si 04 grains are a viable wind scenario for M-type AGB stars, given the 
success of the current models in reproducing observed mass-loss rates, wind velocities, and photometry. Both synthetic and observed 
photometry suggest that the dusty envelopes of M-type AGB stars are quite transparent at visual and near-IR wavelengths, otherwise 
the variations in visual flux would not be dominated by molecular features. 

Key words. Stars: late-type Stars: AGB and post-AGB Stars: atmospheres Stars: mass-loss Stars: winds, outflows, circumstellar 
matter, dust 


1. Introduction 


It is generally assumed that the slow but massive outflows of gas 
and dust from cool luminous AGB stars are accelerated by radi¬ 
ation pressure on dust grains. The total momentum of the pho¬ 
tons emitted by these stars easily matches, or even exceeds, the 
typical momentum of the stellar winds, and solid particles with 
the right optical properties can be very efficient in gaining mo¬ 
mentum from stellar photons through absorption and scattering 
processes. 


Observationally, grain materials are usually identified 
through their characteristic lattice modes in the mid-IR, e.g. 
the well-known silicate features at about 10 and 18 pm which 
are due to stretching and bending modes in the SiC >4 tetrahe¬ 
dron. Mid-IR spectra of AGB stars give valuable insights into 
the complex dust chemistry in the circumstellar envelopes (see 
e.g. Dorschner||20T0 ; Molster et al. 2010 for an overview). For 
determining which grains may contribute to wind acceleration, 
however, important constraints come from the visual and near- 


infrared wavelength regions which are cmcial for the energy and 
momentum budget of the atmosphere since the radiative flux 
maximum of the star is at about 1-2 pm. 

This paper is the third in a series dedicated to identifying 
wind-driving dust species in M-type AGB stars by using a com¬ 
bination of different dynamical models and observational data. 
In the first paper (Bladh & Hofner] |20 1 2| > we focused on dy¬ 
namical constraints for material properties, using a simple pa- 
rameterised description of the dust component in frequency- 
dependent radiation-hydrodynamical models for pulsating atmo¬ 
spheres and winds. Since dust temperatures are strongly affected 
by the wavelength-dependence of the grain opacities (causing 
greenhouse or inverse greenhouse effects), we found that many 
dust species cannot condense sufficiently close to the star to trig¬ 
ger or accelerate a wind. In particular, Fe-bearing silicates suffer 
from a severe greenhouse effect as a result of the steep slope of 
the absorption coefficient at near-IR wavelengths, moving their 
condensation zone much further away from the star than for 
Fe-free magnesium silicates, which was also demonstrated by 
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|Woitke| ( f2006j ). This effect, by definition, does not appear in wind 
models with grey radiative transfer, leading to a severe underes¬ 
timation of the condensation distance for Fe-bearing silicates. 

In the second paper ( (Bladh et al.[ |2013| ) we presented syn¬ 
thetic spectra and visual and near-IR photometry resulting from 
these parameterised models (set P) and from more detailed wind 
models driven by photon scattering on Mg 2 SiC >4 grains (set D), 
first presented in Hofner ( [2008b[ >. Comparing the synthetic pho¬ 
tometry with observations of M-type AGB stars we found that 
the latter models (set D) give a good agreement. In particular, 
they can reproduce large variations of the (V - K ) colour dur¬ 
ing the pulsation cycle due to abundance variations of TiO. This 
is a strong indication that the circumstellar envelopes of M-type 
AGB stars, and, consequently, the wind-driving dust grains are 
rather transparent, excluding true absorption by dust as the main 
source of momentum. Results from the parameterised models 
(set P) support this conclusion, which makes photon scattering 
on Fe-ffee silicates a prime candidate for the wind driving mech¬ 
anism. In order for the grains to be efficient at scattering this 
scenario requires particles of sizes in the range of 0.1 - 1 pm. 
Recently Norris et al. (2012) claimed the detection of grains 
with sizes of ~ 0.3 pm in the close circumstellar environment 
of 3 AGB stars, using multi-wavelength aperture-masking po- 
larimetric interferometry. 

Given the success in reproducing dynamic and photometric 
properties for the small sample of wind models driven by photon 
scattering on Mg 2 Si 04 grains, we here present the first extensive 
set of time-dependent wind models for M-type AGB stars based 
on this mechanism. The aim of this study is to demonstrate that 
outflows can be produced not only for a few selected cases, but 
also for models with a wide range of stellar parameters. We also 
investigate how the observable properties of these models are 
affected by various parameters, and we compare our results with 
observations, in particular dynamic properties, visual and near- 
IR photometry, and mid-IR spectra. 

The paper is organised as follows: in Sect.[2]we introduce the 
basic physical assumptions and input parameters for the dynam¬ 
ical models of M-type AGB stars. In Secs. m we evaluate the 
wind properties, the visual and near-IR photometry and the mid- 
IR spectra produced by these models and compare with observed 
values. In Sect. [6] we explore constraints on the position of the 
wind acceleration zone and in Sect. [7] we provide a summary of 
our conclusions. 


nanometre range, followed the much slower growth of the grains 
through condensation of material from the gas phase, on time 
scales comparable to atmosphere and wind dynamics. Currently 
there is no well-established nucleation theory for oxygen-rich 
atmospheres, describing the formation of the seed particles. 
Fundamental questions, e.g. their composition (which may differ 
from the bulk of material building up the grains) are still a mat¬ 
ter of debate. Our modelling of dust formation follows Gail & 
Sedlmayr ( |1999| , assuming that the growth rate of Mg 2 Si 04 is 
determined by the addition of SiO molecules to the grain surface, 
and using a parameter lu/nH that sets the number of seed parti¬ 
cles per hydrogen atorrPl These seed particles will start to grow 
when the thermodynamic conditions are favourable. The equa¬ 
tion describing the time-dependent growth of Mg 2 Si 04 grains 
is discussed in |Bladh et al.| ( (20131 1. Since we assume that the 
formation of seed particles precedes grain growth (i.e. no new 
seed nuclei are formed) grain size will be uniform in a given 
mass layer. Because of changing conditions in the atmosphere 
and wind, however, grain sizes will differ from layer to layer and 
change with time. 

Each dynamic model is represented by a sequence of snap¬ 
shots of the radial structures, covering typically hundreds of pul¬ 
sation periods. From this long sequence we select a series of 
snapshots, equidistant in phase, during three consecutive pulsa¬ 
tion periods. For these snapshots we produce spectra and pho¬ 
tometry in an a posteriori radiative transfer calculation. The syn¬ 
thetic spectra are computed with opacities from the COMA code 
(Aringer 2000j l and mean synthetic photometry is calculated by 
phase-averaging over the pulsation cycle. 

The current study includes 139 solar-mass models of M-type 
AGB stars, with effective temperatures ranging from 2600 to 
3200K. Three different luminosities have been used: 5000 L 0 , 
7000 L 0 and 10000 L G . These parameters were chosen to repre¬ 
sent common M-type AGB stars, according to stellar evolution 
models. The pulsation period controlling the variations at the in¬ 
ner boundary is derived from the period-luminosity relation pre¬ 
sented in Feast et al. ( | 1989( 1 and the piston velocity amplitude 
Au p ranges between 2-4 km/s, depending on luminosity. This 
results in shock amplitudes of about 15-20 km/s in the inner at¬ 
mosphere. Table|I]shows the combination of input parameters of 
the models, sorted by luminosity. For each fixed luminosity, the 
effective temperature, piston velocity amplitude and seed parti¬ 
cle abundance are varied according to the values listed. 


2. Modelling method and parameters 

The models presented in this paper are similar to model set D 
Bladh et al. (2013)1, corresponding to the models first pre¬ 


in 


sented in|Hofner (|2008b[). We only give a short summary of their 


main ingredients here. For more detailed information concerning 
the dynamical models and the a posteriori radiative transfer, see 
previous artic les (|Hofner et al.| |2003( |Hofner| |2008b( |Bladh &| 
|Hofner([2012HBladh et al.||2013) and references therein. 

The variable structures of the atmospheres and winds are 
produced by simultaneously solving the equations of hydro¬ 
dynamics, frequency-dependent radiative transfer and time- 
dependent grain growth, assuming spherical symmetry. The stel¬ 
lar pulsations are simulated by adopting sinusoidal variations of 
velocity and luminosity at the inner boundary of the atmosphere. 
When calculating the radiative acceleration of Mg 2 Si 04 parti¬ 
cles, grain-size dependent opacities are used, taking into account 
both scattering and absorption effects. 

Dust formation is considered as a two step process, start¬ 
ing with the rapid formation of tiny seed nuclei with sizes in the 


3. Wind properties 

A schematic overview of the dynamical properties of the avail¬ 
able models is shown in FigJT The range in effective temperature 
and luminosity is listed on the x-axis and y-axis, respectively. 
Each subset consists of twelve boxes and is organised such that 
the piston velocity is increasing upwards and the seed particle 
abundance is increasing towards the right. The red boxes repre¬ 
sent models that develop a stellar wind, the blue boxes represent 
models where no wind forms and the grey boxes indicate com¬ 
binations of parameters not tested or where the models fail for 
numerical reasons. It is clear from this plot that the dynamical 
models can produce outflows for a wide range of stellar parame¬ 
ters, although it is generally more difficult to drive winds for high 


1 In contrast, nucleation in carbon-rich stars is often described by 
classical nucleation theory, assuming that the seed particle and micro¬ 
scopic grains consist of the same material, i.e. amorphous carbon. The 
values used here for the parameter n gr /n H are comparable to values 
found in carbon star models. 
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Table 1 . The combinations of input parameters for the models of M-type AGB stars, sorted according to luminosity. The columns 
list stellar mass M*, stellar luminosity L*, effective temperature 7’* of the star, pulsation period P, piston velocity amplitude Am p 
and the assumed seed particle abundance n gr /ttH- The pulsation period is derived from the period-luminosity relation presented in 
|Feast et al.| ( | 1989| >. 


Model series 

M* 

[47 0 ] 

u 

[G/l 

r* 

[K] 

P(U) 

[d] 

A Up 
[km/s] 

log ngjiia 
[a, b, c, d] 

L50 

1 

5000 

2600-3000 

310 

3.0, 4.0 

-16.0, -15.5, -15.0,-14.5 

L70 

1 

7000 

2600-3100 

395 

3.0, 4.0 

-16.0,-15.5,-15.0,-14.5 

L10 

1 

10000 

2600-3200 

525 

2.0, 3.0 

-16.0,-15.5,-15.0,-14.5 



Fig. 1 . Dynamic behaviour of the models as a function of in¬ 
put parameters. The red rectangles indicate models with a stellar 
wind, the blue rectangles indicate models where no wind forms, 
and the grey rectangles indicate combinations of parameters not 
tested or where the models fail for numerical reason. For each 
combination of luminosity and effective temperature the seed 
particle abundance and piston velocity are varied as indicated by 
the inset box. See Table[l]for the different values of seed particle 
abundance n SI /nn (a-d). 


effective temperatures and low luminosities. Judging from these 
trends, models with lower effective temperature and/or higher lu¬ 
minosity than in this current sample can be expected to produce 

outflows^ 

While the effects of certain parameters are rather straight¬ 
forward (L ,, r„), or have been discussed in other papers (Aw p ), 
the seed particle abundance needs some further discussion. The 
range of values for this parameter is limited on both ends since 
photon scattering as a wind-driving mechanism requires grains 
in a certain size range (see Hofner 2008a for a discussion). 
Fewer seeds make the growth of individual grains more effi¬ 
cient (less competition for material), but may actually lead to 
a smaller collective opacity if the grains grow beyond the op¬ 
timal size range. More seeds lead to stronger competition for 
condensable material, resulting in smaller grains. This effect is 


2 The current set of models for M-type AGB stars does not lead to any 
wind models with episodic outflows, whereas such models are found in 
the grid for C-type AGB stars of Mattsson et al.| ( [2010| >, as discussed 
by Eriksson et al. (20141. A possible reason for the different wind be¬ 
haviour is that dust particles with a large absorption cross-section in the 
visual and near-IR (e.g. amorphous carbon) cause back-warming effects 
in the atmosphere. This back-warming can influence the dust formation 
in the layers below and consequently the dynamics. If the wind-driving 
dust species is very transparent in the visual and near-IR (e.g. mag¬ 
nesium silicates) and acting predominately through photon scattering, 
then this feedback through back-warming will not occur. 


demonstrated in Fig. [2] showing the grain radius averaged over 
pulsation cycle in the outermost layers of the dynamical mod¬ 
els that produce a wind, sorted according to seed particle abun¬ 
dance. The degrees of condensed silicon (Si) in the models that 
produce outflows are shown in Fig. [3] with typical values rang¬ 
ing between 10-25%, corresponding to dust-to-gas mass ratios 
of about 10 2 . This is in reasonable agreement with recent re¬ 
sults by Khouri et al. (2014) for W Hya, who find that about a 
one third of Si is condensed into grains. According to our models 
the degree of condensed magnesium will be about twice as much 
as for Si since Mg 2 SiC >4 has two magnesium atoms per silicon 
atom and the element abundances of Si and Mg are comparable. 

A comparison of wind velocities and mass-loss rates of 
the models with corresponding empirical data for M-type AGB 
stars, derived from observations of several CO-lines, is shown in 
Fig. 0 The agreement with observations is good, not the least 
considering that this set of models does not correspond to a stel¬ 
lar population, but rather a grid of input parameters where not 
all combinations of parameters are equally likely. Both high and 
low mass-loss rates are reproduced, and higher mass-loss rates 
could probably be reached by increasing the stellar luminosity 
of the wind models further. 

The four panels of Fig. [4] show the same data but the mod¬ 
els are colour-coded in different ways to illustrate the effects of 
different input parameters. The top left panel of Fig. [4] shows 
the dynamical properties of the wind models, colour-coded ac¬ 
cording to luminosity. The models with different values for the 
luminosity form bands where higher luminosity correlates with 
higher mass-loss rates. If we instead plot the dynamical proper¬ 
ties colour-coded by seed particle abundance (lower left panel), 
the models form bands correlating more with wind velocity. In 
right panels of Fig.[4]we also plot the models sorted according to 
effective temperatures and piston velocities, but we see no dis¬ 
tinct trends in wind properties with respect to these parameters. 


4. Visual and near-IR photometry 

Dust-driven winds can be accelerated by both absorption and 
scattering of stellar photons on dust grains, but scattering is only 
significant if the grains grow to adequate sizes. If true absorption 
is the dominant source of the radiative acceleration then sig¬ 
nificant circumstellar reddening will occur; the dusty envelope 
will veil the molecular features in the visual wavelength region 
and emit a significant infrared excess. If scattering is the dom¬ 
inant source of the momentum there will be much less circum¬ 
stellar reddening; changes in the molecular features during the 
pulsation cycle will be visible at visual wavelengths and there 
will not be much infrared excess, even though dust is present. 
Since the visual and near-IR spectra are strongly sensitive to the 
wind-driving mechanism, i.e. the optical properties of the wind¬ 
driving dust species, we use the colours ( J-K ) and (V-K) when 
comparing model results with observational data. The sources of 
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Fig. 4. Observed mass-loss rates vs. wind velocities of M-type AGB stars (Olofsson et al. 


2002 


Gonzalez Delgado et al. 


2003 


plus 

signs) and the corresponding properties of all models which develop a wind (indicated in red in Fig. |T]». All panels show the same 
data but with different colour coding according to stellar luminosity (upper left), effective temperature (upper right), seed particle 
abundance (lower left) and piston velocity amplitude (lower right). The stellar mass is one solar mass for all models. 


Table 2. Input parameters for the dynamical models resulting in 
the photometric variations shown in Fig. [5] 
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the observational data sets used here are described in detail in 
|Bladh et al.l ( |20T3) . 

AGB stars are variable stars with photometric properties that 
change with phase. Realistic dynamical models should be able 
to reproduce this photometric variation during the pulsation cy¬ 
cle. In Fig. 5] we plot the synthetic colour (J - K ) vs. (V - K) 
for a few se ected models during a pulsation cycle, together with 
phase-dependent colours constructed from sine fits of observed 
light curves (see Bladh et al. 2013 for more details). These 
photometric variations form loops in the colour-colour diagram. 


with different positions, shapes and tilts. For M-type AGB stars, 
both the synthetic and observed colour loops are characterised 
by large variations in (V - K ) and small variations in (J - K ), as 
can be seen in Fig. [5] The colour (V - K) reaches its peak value 
during luminosity minimum, followed by a trough during lumi¬ 
nosity maximum. A closer examination of the synthetic spectra 
reveals that the large variations in (V - K ) during a pulsation cy¬ 
cle are caused by a change in the abundance of molecules (TiO, 
see Fig. [6]), and not by changes in the amount of dust present 
(Bladh et al. |2013) l. This is a strong indication that the wind- 
driving dust species in M-type AGB stars are quite transparent 
in the visual and near-IR, otherwise these variations would not 
be dominated by molecular features. 


Looking at Fig. [5] it is clear that phase-averaged colours will 
be affected by both the position of the loops in the colour-colour 
diagram and the variations in colour during a pulsation cycle. 
Fig-0 shows phase-averaged synthetic colours (filled squares), 
as well as observed colours for both M-type (grey and red plus 
signs) and C-type AGB stars (green plus signs). The synthetic 
colours are characterised by blue values in (J - K ), especially 
compared to values of C-type AGB stars, with small differences 
between the individual dynamical models. There is a strong cor¬ 
relation between the effective temperatures of the models and the 
values of (J - K), as can be seen in the top right panel of Fig. 0 
This is probably because the temperature of the star has a strong 
impact on the water features (see Fig.|6|. 
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Fig. 2. Grain sizes (averaged over pulsation cycle) at the outer 
boundary of the models that produce a stellar wind, sorted ac¬ 
cording to different seed particle abundance n gI /nu. 



fc 


Fig. 3. The fraction of Si atoms condensed into grains, averaged 
over pulsation cycle at the outer boundary of the models that 
produce a stellar wind. 


The phase-averaged synthetic (V - K ) colours show large 
differences between individual models. The values of the mean 
(V - K) colours are mostly determined by the range of varia¬ 
tion during the pulsation cycle, i.e. how far the loops extend to 
the blue side of the diagram. As mentioned above, the variation 
in visual fluxes is determined by changes in TiO, which in turn 
reflects changes in temperature during the pulsation cycle. The 
stronger the temperature variations, the wider the loops, result¬ 
ing in bluer average (V - K) colours. 


The conclusion we can draw from Fig. [7] is that the mean 
synthetic colours reproduce well the observed colours from field 
M-type AGB stars (Mendoza, \961) and the bulk of the observed 
values from Bulge miras (Groenewegen & Blommaert 2005) >, 
especially the large spread in (V - K). The observed colours 
from Bulge miras show a much larger spread in (J - K) than 
the synthetic colours, but this could partly be due to metallicity 
effects and the fact that we assume solar abundances of C and O 
in the models. A higher C/O-ratio (closer to 1 than the 0.48 that 
we assume) should give redder values of (V - K ) and (J - K ), 
judging from hydrostatic model atmospheres. 



Fig. 5. Photometric variations during a pulsation cycle in (J - K ) 
vs. (V-K). The black loops show photometric variations derived 
from sine fits of observed light-curves for a set of M-type miras 
(R Car, R Hya, R Oct, R Vir, RR Sco, T Col and T Hor) and 
the coloured loops show synthetic photometric variation from a 
few selected models (see Table [2}. For the observed miras we 
adopted photometric data in the visual from Eggen ( 1975| l and 
Mendoza] ( 1967|> and complemented those with the near-IR data 


published by Whitelock et al. 


calculated from sine fits of the 


2000)1. The synthetic loops are 


ight-curves in the same way as 


the observational data. See [Bladh et al.| ( |2013| ) for more details 
concerning the sinefit of observed light curves. The bottom panel 
shows the same data as the top panel, but zoomed in. The obser¬ 
vations in grey (plus signs) are presented in Fig. [7] 



X [pm] 

Fig. 6. Spectral energy distribution as a function of wavelength 
for a model with input parameters M = 1M 0 , L = 5000 L 0 , 
T * = 2800K, u p = 4km/s, and logn gr /nH = 14.5. The original 
spectrum is shown in black and the coloured curves show spectrg 
calculated without TiO opacities (red) and FEO opacities (blue), 
indicating which photometric bands will be affected by changes 
in the corresponding features. 
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Fig. 7. Observed and mean synthetic colours (J - K ) vs. (V - K ). The observational data (plus signs) is compiled from different 
sources: Galactic Bulge miras ( |Groenewegen & B lommaert 2003] grey), field M-type LPVs (Mendoza 1967, red), and C-rich giants 
(Bergeat et al. 2001 [ green). The model data (filled squares) is from models which develop a wind, as indicated in red in Fig. [T| All 
panels show the same data but with different colour coding according to stellar luminosity (upper left), effective temperature (upper 
right), seed particle abundance (lower left) and piston velocity (lower right). The stellar mass is one solar mass for all models. Most 
of the observational data are single epoch measurements, whereas the synthetic colours are means over the pulsation cycle. 


5. Mid-IR spectra 

An issue with the current wind models is that, although they 
produce dynamical properties and photometry in the visual and 
near-IR in agreement with observations (see Secs.[3]an d 0 , they 
do not show the prominent silicate features at 10 pm and 18 pm 
that are observed in many oxygen-rich AGB stars. These sili¬ 
cate features have been identified as belonging to magnesium- 
iron silicates, such as olivine ([Mg,Fe] 2 Si 04 ) and pyroxene 
([Mg,Fe]Si 03 ), and are caused by stretching and bending res¬ 
onances in the Si 04 -tetrahedron. 

The optical properties of Fe-free silicates offer an explana¬ 
tion as to why the synthetic spectra do not show these promi¬ 
nent features. In the wavelength region where AGB stars emit 
most of their flux the absorption cross-section of magnesium sil¬ 
icates is very low, in particular compared to the spectral fea¬ 
tures at 10 pm and 18 pm, so the heating of the grains by 
stellar light is moderately efficient. Fig. [ 8 ] shows the absorbed 
and emitted light for Mg 2 Si 04 grains (black curves), assuming 
a geometrically diluted radiative flux from a hydrostatic atmo¬ 
sphere (with stellar parameters typical of an M-type AGB star 
M = 1M 0 , L = 5000L g and T, = 2800 K) and a dust tempera¬ 
ture of 7 d = 1000K. The energy balance is set mostly in the 
mid-IR wavelength region and the grains cool efficiently. This 
is the reason why such grains can form close to the stellar sur¬ 


face, but it also leads to a rapidly falling dust temperature in the 
wind. 

In order to produce silicate features, the grain temperature 
needs to fall less rapidly with distance from the star than for 
pure Mg 2 Si 04 grains. Fig. [ 8 ] indicates that the inclusion of Fe 
in the silicates will help to achieve this; for MgFeSi 04 grains 
(red curves) light is mostly absorbed in the near-IR wavelength 
region and then emitted in mid-IR region. However, because of 
the increased near-IR absorption Fe-bearing silicates can only 
become thermally stable further away from the star. Simple 
estimates and results from dynamical models with frequency- 
dependent radiative transfer predict that silicates with significant 
Fe content (e.g. MgFeSKIL) become stable at 5 - 10 R* (jWoitke] 
[2006] |Hofner| |2008a[ |Bladh & Hdfner|[20T2| . |Woitke1 ( |2006| l in- 
vestigated the growth of inhomogeneous dust grains, composed 
of a mixture of Mg 2 Si 04 , Si 02 , ALO 3 , Ti 02 , and solid Fe in a 
stellar outflow. He found that once Fe inclusion is possible, it 
will control the temperature of the inhomogeneous dust grains 
and act as a thermostat, keeping the grain temperature just be¬ 
low the sublimation limit. The amount of included Fe is thereby 
adjusted in a self-regulating process. 

To test this scenario we use a snapshot of the radial structure 
from a dynamical model where the wind is driven by photon 
scattering on Mg 2 Si 04 grains. In the a posteriori radiative trans- 
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Fig. 8. Absorbed and emitted flux, Q ahs (d)F ( and <2abs(d-)^.u 
for two different grain materials, Mg 2 SiC >4 (black) and 
MgFeSiOa (red). The absorbed light is calculated using a ge¬ 
ometrically diluted radiation field from a hydrostatic atmo¬ 
sphere (at 2R* for MgoSiC^ and at 10 R* for MgFeSKTt) 
with stellar parameters typical of an M-type AGB star 
(M = 1M 0 , L = 5000 L 0 and T t = 2800 K). The emitted light is 
calculated using a Planckian source function at a grain temper¬ 
ature of I'd = 1000 K. The optical data for Mg 2 SiC >4 is from 
Jager et al. (2003[i and the optical data for MgFeSKTt is from 


Dorschner et al.|(ji995|l. We note that the local radiation field in 


the wind-acceleration zone will be more pronounced in the mid- 
IR region than what is shown here because of the addition of the 
emitted light from the dust component. 


fer calculations we set a lower limit on the grain temperature at 
800 K or 500 K, as if the grains contained impurities that would 
keep them warm. The top panel of Fig.[9]shows the original tem¬ 
perature structure for the Mg 2 SiC >4 particles, as well as the cases 
where a minimum dust temperature of 800 K and 500 K is set. 
The resulting spectra are shown in the top panel of Fig. [TO] It is 
clear from the prominent silicate features produced in the latter 
cases that the missing spectral features are not caused by a low 
abundance of silicate dust in the dynamical models, but rather a 
too low grain temperature. Keeping the grains warm results in 
silicate emission features. The warmer the grains are, the more 
prominent the features. 

As mentioned above, a way to increase the grain tempera¬ 
ture is to include impurities of Fe in the Mg 2 SiC >4 grains. This 
can be simulated by considering dust particles with a core of 
Mg 2 SiC >4 and a thin mantle of MgFeSKTt, as if a small frac¬ 
tion of Fe had condensed onto the surface of pure magnesium 
silicate grains^] We determine the temperature of these grains 
through radiative equilibrium, assuming that the addition of a 
thin mantle MgFeSiC >4 will not significantly change the struc¬ 
ture of the dynamical model. In the a posteriori spectral calcu¬ 
lation we replace the temperature structure and the optical prop¬ 
erties of the Mg 2 SiC >4 grains with the corresponding properties 
for the core-mantle grains from the point where they are ther¬ 
mally stable. We assume this to be at the distance from the star 
where the temperature of the core-mantle grains drops below 

1 It is possible to introduce this modification in the a posteriori ra¬ 
diative transfer since the inclusion of a small amount of other materials 
into these particles should not alter the overall structure of the mod¬ 
els significantly; transparent materials will not have enough radiative 
cross-section to alter the energy balance and opaque materials are only 
thermally stable further out in the wind. 
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Fig. 9. Grain temperature as a function of distance from the star 
for a model with stellar parameters M = 1 M G , L = 5000 L G and 
T t = 2800 K during luminosity minimum. Top panel: The origi¬ 
nal temperature structure of the Mg 2 Si 04 grains (black) and the 
temperature structure when setting a minimum temperature of 
800 K (red) and 500 K (blue). Lower panels: The original tem¬ 
perature structure of the Mg 2 Si 04 grains (black) and the tem¬ 
perature structure of dust particles with a core of Mg 2 Si 04 and a 
mantle of MgFeSi 04 when it is thermally stable (T < 1000 K), 
as indicated by the vertical grey line. The green, red and blue 
lines show the temperature structure for a mantle size that is 1%, 
5% and 10% of the grain radius, respectively. 



1000K (approximately the temperature where magnesium-iron 
silicates will condense). The resulting temperature structure for 
three different mantle sizes (1%, 5% and 10% of the radius cor¬ 
responding to 3%, 15% and 30% of the volume) are shown in 
Fig . [9] and the corresponding spectra are shown in the bottom 
panel of Fig. 10 We note that dust particles with a thin mantle 
of MgFeSi 04 may be thermally stable quite close to the stellar 
surface compared to pure MgFeSi 04 grains. For the minimum 
luminosity phase of this particular model the core-mantle grains 
become thermally stable at ~ 2.5 R„, ~ 4.5 R, and ~ 5.5 R*, re¬ 
spectively, if the mantle is 1%, 5% or 10% of the grain radius. 
The overall temperature structure and, consequently, the radial 
distances at which dust particles are thermally stable, will vary 
during the pulsation cycle, and depend on the stellar parameters. 


To summarise the results from the second test, it seems that 
a thin mantle of Fe-bearing silicates, in this case consisting of 
MgFeSi 04 , can be thermally stable quite close to the stellar 
surface. This thin mantle of MgFeSKT; corresponds to a small 
Fe/Mg-ratio but is enough to heat the dust particles and pro¬ 
duce silicate features without changing the spectra in the visual 
and near-IR significantly. For illustration. Fig [TT] shows the ISO 
spectra of two M-type AGB stars, R Aqr and R Cas, in addition 
to the original model spectrum and the model spectra with a thin 
mantle of MgFeSi 04 on top of the Mg 2 Si 04 grains. 
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Fig. 10. Spectral energy distributions as a function of wave¬ 
length resulting from the dust temperatures indicated in Fig. [9 
The black curve shows the original spectrum for a dynamica 
model with pure Mg 2 SiC >4 grains (M = 1 Mg, L = 5000L G and 
r* = 2800 K). Top panel: The spectra produced when setting a 
minimum grain temperature of 800 K (blue) and 500 K (red) for 
the pure Mg 2 SiC >4 grains, i.e. using the temperature structure 
in the top panel of Fig. [9] Bottom panel: The spectra produced 
when using a core of Mg 2 SiC >4 and a mantle of MgFeSiC> 4 . The 
green, red and blue lines show the spectra for a mantle size of 
1%, 5% and 10% of the radius, respectively. The corresponding 
temperature structures are shown in the lower panels of Fig. [9] 


6. Grain temperature and the onset of outflow 

The grain temperature in our dynamical models is assumed to 
be determined by radiative equilibrium, but the low absorption 
cross-section of Mg 2 SiC >4 in the wavelength region where AGB 
stars emit most of their flux raises the question if this assump¬ 
tion is realistic. Gauger et al. ( |i990| have estimated the relative 
importance of different processes (radiative, collisional, chem¬ 
ical) for changing the internal energy of dust particles consist¬ 
ing of amorphous carbon or astronomical silicates. They found 
that radiative processes dominate by many orders of magnitude 
over collisional processes in both cases. However, these mate¬ 
rials have high absorption cross-sections in the visual and near- 
IR wavelength regions, leading to efficient radiative heating. In 
Appendix [A] we give estimates for the radiative and collisional 
heating rates of Mg 2 Si 04 grains: while radiative heating is much 
less efficient than for the materials tested by Gauger et al. ( jl990| l, 
it still clearly dominates over heating by collisions with the 
warmer gas. 



10 
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Fig. 11. Spectral energy distribution as a function of wavelength. 
The red curves shows the observed ISO/SWS spectra of R Aqr 
and R Cas ( |Kraemer et al.| |2002[ |Sloan et al.| |2003| l. The black 
solid curve shows the spectra from a dynamical model using 
pure Mg 2 Si 04 grains and the grey curves show the spectra when 
assuming a core of Mg 2 Si 04 and a mantle of MgFeSi 04 in the a 
posteriori spectral calculations. The dashed, dot-dashed and dot¬ 
ted lines show the spectra setting the mantle size to 1%, 5% and 
10% of the grain radius, respectively. 


Under radiative equilibrium conditions, grains of this mate¬ 
rial tend to cool rapidly in the wind (because of their optical 
properties, as discussed in Sect. [5j, with grain temperatures gen¬ 
erally a few hundred degrees below the gas temperature. Setting 
the grain temperature equal to the warmer gas temperature, in¬ 
stead of assuming radiative equilibrium, results in shifting the 
position of the dust formation zone slightly outwards. This pro¬ 
vides us with an opportunity to explore how the radial position 
of the wind acceleration zone affects dynamic and photometric 
properties. Fig. 12 shows the position of the onset of outflow, 
averaged over a pulsation cycle, for the two types of models. 
The red area shows the results from models with a grain tem¬ 
perature determined by radiative equilibrium, indicating a typi¬ 
cal distance of the onset of outflows of about 2 -3 R„, whereas 
the blue area shows the results from the models where the grain 
temperature is set equal to the gas temperature, indicating typical 
distances of about 3-6 R* (the purple area shows the overlapping 
region). 

This shift in the radial position of the wind acceleration 
zone affects the observable properties of the dynamical models. 
Fig.[l3]shows the mass-loss rates and wind velocities of the mod¬ 
els with a dust temperature based on radiative equilibrium (red 
squares) and the models where we have assumed that the dust 
temperature is equal to the gas temperature (blue squares). Both 
fit observations reasonably well, but the blue models (where the 
wind acceleration zone is shifted outwards to lower densities) 
do not reproduce the observations with high values of mass-loss 
rate and wind velocity, while the stars with low mass loss and 
slow wind velocities seem to be better reproduced. 

In contrast to the dynamical properties, photometry provides 
strong constraints on the position of the wind acceleration zone. 
The top panel of Fig. 14 shows the average colours (J - K ) and 
(V - K ) during a pulsation cycle for the two types of models. 
The photometry of the models where the grain temperature is 
set equal to the gas temperature (blue squares) is very similar to 
the photometry for pulsating models without wind (see Fig. 4 in 
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Fig. 12. Position of the onset of outflow for all models devel¬ 
oping a wind, averaged over the pulsation cycle. The red and 
blue areas represent models where the grain temperature is de¬ 
termined by radiative equilibrium and models where the grain 
temperature is set equal to the gas temperature, respectively. The 
purple area is where the two types of models overlap. 


|Bladh et aL| ( |20i3| >), whereas the photometry of the models based 
on radiative equilibrium (red squares) are bluer in (V-K). An in- 
depth comparison is given in the lower panel of Fig.[l4j showing 
the photometric variations during a pulsation cycle. The red and 
blue loops are from two dynamical models where the dust tem¬ 
perature is set in different ways but which otherwise have the 
same model parameters. The black loops show observed pho¬ 
tometric variations from a set of very regular mira stars (for a 


more detailed description of the observational data see Bladh 
et al. (2013|>). The large variations in (V - K ) seen in the ob¬ 


served photometry are not reproduced in the models with a grain 
temperature set equal to the gas temperature. 

It should be mentioned here again that the effects on the pho¬ 
tometric colours are not due to dust opacities. Instead the visual 
and near-IR photometry of the models gives direct diagnostics 
of the gas, i.e. the structure and dynamics of the atmosphere. As 
mentioned in Sect|4] the large variations observed in (V - K) are 
a consequence of a change in the abundance of TiO during the 
pulsation cycle. The position of the onset of outflow affects the 
overall density and temperature structure of the gas and, con¬ 
sequently, the molecular abundances, by radially expanding the 
atmosphere outwards. If the wind acceleration zone is shifted 
to distances where this effect no longer influences the structure 
of the inner atmosphere, then the resulting photometry will re¬ 
semble that of a pulsating atmosphere without a stellar wind. It 
seems that the observed variations in (V - K) do not only reveal 
that the wind-driving dust species in M-type AGB stars have to 
be quite transparent in the visual and near-IR, but they also pro¬ 
vide an upper limit for the distance of the onset of outflow. We 
note that this upper limit probably also holds if a different mech¬ 
anism than radiation pressure on dust accelerates the wind, since 
the diagnostic presented here is based on gas properties. 


7. Summary and conclusions 

In this paper we present the first extensive set of time-dependent 
models for atmospheres and winds of M-type AGB stars. The 
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Fig. 13. Observed mass-loss rates vs. wind velocities of M-type 
AGB stars 
plus signs) 

models where the dust temperature is determined by radiative 
equilibrium (red squares) or set equal to the gas temperature 
(blue squares). The filled squares mark the dynamical models 
where a more detailed analysis is conducted. 
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and the corresponding properties for the dynamica 




Fig. 14. Observed and synthetic visual and near-IR colours. 
Top panel: Synthetic phase-averaged colours from wind mod¬ 
els where the dust temperature is determined by radiative pro¬ 
cesses (red squares) and the gas temperature (blue squares). 
Bottom panel: Photometric variations derived from sine fits of 
observed light-curves for a sample set of M-type miras. The syn¬ 
thetic photometric variations for the models marked with filled 
squares (M = 1 Mq, L = 5000 L G , T* = 2800 K, w p = 4 krn/s 
and log « gr /«H = -15) in the top panel are also plotted (solid 
red and blue loops). The colours are calculated from sine fits of 
light-curves, same as for the observational data. See [Bladh et a+j| 
( 2013) > for more details concerning the sine-fit of observed light 
curves. For details on the observational data see Fig. [5] 
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stellar winds in these models are driven by photon scattering on 
Mg 2 SiC >4 dust. Grains of this material can form close enough 
to the stellar surface, they consist of abundant materials and, if 
they grow to sizes comparable to the wavelength of the stellar 
flux maximum, can provide sufficient radiative acceleration to 
drive outflows. 

The current study includes 139 solar-mass models, using 
three different luminosities (5000 L 0 , 7000 L 0 and 10000 L G ) 
and effective temperatures ranging from 2600 K to 3200 K (see 
Table [I}. The input parameters describing stellar pulsations and 
seed particle abundance are also varied for each combination 
of stellar parameters, resulting in different atmospheric struc¬ 
tures and wind properties. The models produce outflows for a 
wide range of stellar parameters (see Fig. [TJ, as well as mass- 
loss rates, wind velocities, and visual and near-IR photometry 
in good agreement with observations (see Fig. [4] and Fig. |7J. 
This demonstrates that winds driven by photon scattering on 
Mg 2 SiC >4 grains is a viable scenario for the mass loss of M-type 
AGB stars. However, the current models do not show the char¬ 
acteristic silicate features at 10 and 18 fim. This is not caused 
by a low abundance of silicate dust in the dynamical models, 
but rather a too rapidly falling temperature of Mg 2 Si 04 grains. 
Adding a thin mantle of MgFeSi 04 to these grains, at distances 
where such a mantle can be thermally stable, increases the grain 
temperature and results in spectra with pronounced silicate fea¬ 
tures. 


The visual and near-IR spectra provide constraints on the 
driving mechanism of the outflows. Stellar winds driven by pho¬ 
ton scattering will produce very different spectra than winds 
driven by true absorption, as there is much less circumstellar 
reddening by the dust component in the former case. This is con¬ 
firmed by observed and synthetic colours of M-type AGB stars; 
the visual and near-IR fluxes are characterised by large varia¬ 
tions in (V - K) and small variations in (/ - K) during the pul¬ 
sation cycle (see Fig. [5]>. The variations in (V - K ) are caused by 
changes in the molecular abundances (mostly TiO) and not by 
changes in the dust opacities (see Fig. [6|and ]Bladh etaTl[20T3] >. 
From this we can infer that the envelopes of M-type AGB stars, 
and, consequently, the dust grains, are quite transparent in the 
visual and near-IR wavelengths, otherwise the spectral changes 
would not be dominated by molecular features. The results re¬ 
garding the mid-IR features mentioned above, however, indicate 
that the grains have to include a small amount of impurities (dirty 
silicates). Preliminary tests demonstrate that including a small 
amount of Fe in magnesium silicate grains is sufficient to solve 
the problem of grain temperatures and that the visual and near- 
IR photometry will not be affected significantly. Future models 
should include a more complete treatment of dust formation. 

The models presented in this paper are an important step to¬ 
wards closing a long standing gap in stellar evolution theory. For 
the first time, it is possible to compute both mass-loss rates and 
the corresponding spectra from the same time-dependent atmo¬ 
sphere and wind models of M-type AGB stars. This is essential 
input for models of stellar evolution and population synthesis. 
However, before the model results can be generally applied in 
these contexts it is necessary to compute a bigger grid with bet¬ 
ter coverage of mass, luminosity, effective temperature and other 
fundamental parameters. 
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Appendix A: Heating of Mg 2 SiC >4 grains 

The very low absorption cross-section of Mg 2 Si 04 grains in the 
wavelength region where AGB stars emit most of their flux raises 
the question if grain heating is dominated by radiation or colli¬ 
sion with the warmer gas. Estimates of the processes that deter¬ 
mine the internal energy of a dust particle, and consequently the 
grain temperature, have previously been done by Gauger et al. 
(1990) for amorphous carbon and astronomical silicates. Both 
materials have high absorption cross-sections in the near-IR and, 
consequently, radiative processes dominate over collisional en¬ 
ergy exchange. The mid-IR optical properties of Mg 2 Si 04 grains 
(relevant for radiative cooling, see Sect. [5]i are comparable to 
those of astronomical silicates. However, the near-IR absorption 
coefficients, which determine the radiative heating, are several 
orders of magnitudes lower than those of Fe-bearing silicates. 
Therefore we need to examine if radiative heating still domi¬ 
nates over heating by collisions with the warmer gas in the case 
of Mg 2 Si 04 grains[j 

We start by estimating the energy exchange occurring when 
a dust particle of size a and temperature collides with a gas 
particle of temperature T s , 


AEcoii - x— (T g - Td), 
47r 


Gauger et al. (19901 considered a third process, i.e. energy change 


due to association or disassociation of one monomer. They found this to 
be negligible compared to collisional processes, and similar reasoning 
holds for the current case. 
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(for details see Gauger et al.] p990| >). If we multiply the energy 
transfer during one collision AE C0 \\ with the number of collisions 
per second r“* u , we can estimate the rate of change in internal 
energy for a dust particle as a result of collisions with the sur¬ 
rounding gas, 

1 2 

Tcoii — AEcoil — ttci v’ih ni[/I AE co ]\. 

Tcoll 

Here v* is the average relative velocity between a dust particle 
and a hydrogen atom which, given the difference in mass be¬ 
tween the two particles, can be approximated with the average 
thermal velocity of the hydrogen atomj^The particle density of 
hydrogen is denoted ;; H and [j is the probability that a collision 
will result in an energy transfer of A_E co u. The term t“| m is a mea¬ 
sure of how many gas particles per second collide with the cross- 
section no 1 of a dust particle and transfer energy. If we assume 
that the particle density of hydrogen is given by na - pi 1.4m H 
we get the following expression for the change in energy due to 
collisions 

Lcoii = TTfl 2 I'th/?7X—-pC^g-Td). 

1 ,4oth 47t 

We set the grain temperature to Ta — 1000 K and the gas temper¬ 
ature to T g = 1500 K, both typical values in the dust formation 
zone. Furthermore, we assume a grain size of a = 10 5 6 * * * * cm and 
set p = 1 (corresponding to an upper limit). The average ther¬ 
mal velocity for this gas temperature is v* ~ 5 km/s. According 
to detailed dynamical models, a typical gas density is about 
p ~ 10 14 g/cm 3 in the dust formation zone rj Consequently, the 
energy gained per second due to collisions of a dust particle with 
the surrounding gas is estimated to 

Tcoll ~ 4 ■ 10- 11 [erg/s]. 

Next we estimate the radiative heating rate of a dust particle. The 
rate of change in internal energy of a grain due to the absorption 
of photons from the ambient radiation field can be expressed as 

/-*oo 

Tabs = 47T7R7 2 I Q(A, Cl) J A dd, 

Jo 

where a is the radius of the dust particle, Q(a , A) is the absorp¬ 
tion efficiency and Jj is the mean intensity of the incident ra¬ 
diation field. In the case of M-type AGB stars photometric ob¬ 
servations indicate that the stellar radiation is not significantly 
reddened by the dust material in the circumstellar envelope. We 
can therefore consider the radiation field as optically thin and 
approximate it by a geometrically diluted stellar radiation field, 
here represented by a Planck function. In the small particle limit 
the optical properties of a dust material Q'(A ) = Q(A, a) I a can be 
approximated with a power-law function, where D | cm ‘ 1 ] and p 
are material-dependent properties, 

—r- 

5 The hydrogen particles in the dust formation zone will probably 
mostly be in the form of H 2 , but we ignore this extra complication in the 
following order of magnitude estimates. This means that the following 
numbers will overestimate the collisional heating somewhat, erring on 
the safe side. 

6 This agrees well with observed mass-loss rates and wind veloci¬ 

ties for M-type AGB stars. Using conservation of mass, p = M/4nR 2 v, 

and assuming a typical wind velocity of 10 km/s and mass-loss rate of 

10~ 6 M e /yr results in a density of about p as 10~ 14 g/cm 3 in the dust 

formation zone. 


Consequently, the rate of change in internal energy due to ab¬ 
sorption of photons for a dust particle in the small particle limit 
is given by 


Tabs = 4 ttV 
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where x = hc/kAT. Here D and p are fitted in the wave¬ 
length region where the star emits most of its flux, which is 
approximately at 1 pm for a star with an effective tempera¬ 
ture T, = 2800 K. Using the optical data for Mg 2 Si 04 given 
by J tiger et al. (2003) we derive material constants D « 3 cm 1 


and p ~ -0.9 for do = 1 pm. Since Mg 2 Si 04 grains condense 
around r « 2 R* the geometric dilution factor is approximately 
W(r) = | [J - -^1 - (/?*/r) 2 ] ss 0.07. Setting a = 10~ 6 cm, and 
T(4 + p)(T(4 + p) ~ 2.6 we estimate the rate of change in inter¬ 
nal energy for a dust particle due to absorption of photons to be 


Tabs -2-10 erg/s. 


It should be noted that this number is probably a lower limit 
since the power law fit of Q' used here tends to underestimate 
the contribution to grain heating from the Mid-IR wavelength re¬ 
gion. There the lower stellar flux may be compensated by higher 
grain opacity around 10pm (see Fig. [8j. 

Given these simple estimates we can safely assume that ra¬ 
diative heating dominates over collisional heating for magne¬ 
sium silicates, even though the absorption cross-section of these 
grains is very low in the wavelength range where M-type AGB 
stars emit the most of their flux. These estimates are for quite 
small grains (0.01 pm) but should hold for even larger grains, 
given the functional dependence of each process on grain size. 
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